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Abstract 

The relation between the Star Formation Rate (SFR) and stellar mass (M*) of galaxies represents a fundamental constraint 
on galaxy formation, and has been studied extensively both in observations and cosmological hydrodynamic simulations. 
Flowever, the observed amplitude of the star formation rate - stellar mass relation has not been successfully reproduced in 
simulations, indicating either that the halo accretion history and baryonic physics are poorly understood/modeled or that 
observations contain biases. In this paper, we examine the evolution of the SFR-M* relation of z ~ 1 - 4 galaxies and 
display the inconsistency between observed relations that are obtained using different techniques. We employ cosmolog¬ 
ical hydrodynamic simulations from various groups which are tuned to reproduce a range of observables and compare 
these with a range of observed SFR-M* relations. We find that numerical results are consistent with observations that use 
Spectral Energy Distribution (SED) techniques to estimate star formation rates, dust corrections and stellar masses. On 
the contrary, simulations are not able to reproduce results that were obtained by combining only UV and IR luminosities 
(UV-l-IR). These imply SFRs at a fixed stellar mass that are larger almost by a factor of 5 than those of SED measure¬ 
ments for z~ 1.5-4. For z < 1.5, the results from simulations, SED fitting techniques and IR-l-UV conversion agree 
well. We find that surveys that preferably select star forming galaxies (e.g. by adopting Lyman-break or blue selection) 
typically predict a larger median/average star formation rate at a fixed stellar mass especially for high mass objects, with 
respect to mass selected samples and hydrodynamic simulations. Furthermore, we find remarkable agreement between 
the numerical results from various authors who have employed different cosmological codes and run simulations with 
different resolutions. This is interesting for two reasons. A) simulations can produce realistic populations of galaxies 
within representative cosmological volumes even at relatively modest resolutions. B) It is likely that current numerical 
codes that rely on similar subgrid multiphase Inter-Stellar Medium (ISM) models and are tuned to reproduce statistical 
properties of galaxies, produce similar results for the SFR-M* relation by construction, regardless of resolution, box 
size and, to some extent, the adopted feedback prescriptions. 
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1 Introduction 

Numerous studies have demonstrated a correlation between 
galaxy star formation rates and stellar masses at redshifts 
z ~ 1-4. These studies have employed a range of different 
observing strategies to sample the galaxy population as com¬ 
pletely as possible (Noeske et al. 2007; Elbaz et al. 2007; 
Daddi et al. 2007; Drory & Alvarez 2008; Kajisawa et al. 
2010; Karim et al. 2011; Bauer et al. 2011; Gilbank et al. 
2011; Reddy et al. 2012; Bouwens et al. 2012; Heinis et al. 
2014; Koyama et al. 2013; Guo et al. 2013; de Barros et al. 
2014; Speagle et al. 2014). Galaxies in different surveys and 
redshifts are usually selected using different techniques and 
wavelengths. Some examples of selecting galaxies in the lit¬ 
erature are: 
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1. Lyman-break selection: in the absence of dust extinc¬ 
tion, star forming galaxies have a flat continuum at rest 
frame ultraviolet wavelengths. However, blueward of 
the Lyman break (rest-frame 912 A), photoelectric ab¬ 
sorption by galactic or extragalactic sources of neutral 
hydrogen sharply cuts the emitted spectrum. High red- 
shift sources can be selected using this spectral break. 
This method has been extensively used to construct 
samples of galaxies at different epochs and study their 
Star Eormation Rate-Stellar Mass (SFR-M*) relation 
(Reddy et al. 2012; Bouwens et al. 2012; Heinis et al. 
2014; de Barros et al. 2014). An important drawback 
of this technique is that it is only capable of selecting 
objects with young ages, high star formation rates 
and low dust content (i.e star forming galaxies). Only 
these objects are able to produce a large amount of UV 
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light which is not absorbed by dust. This makes the 
Lyman-break technique insensitive to galaxies that are 
dust-reddened or contain evolved stellar populations. 

2. Ha selection; in this case, galaxies are selected by their 
emission in the Ha line (rest-frame 6563 A), which is 
correlated with star formation. Various authors have 
constructed Ha selected samples to study the evolution 
of the cosmic star formation rate density and SFR-M* 
relation (e.g. Sobral et al. 2013; Koyama et al. 2013; 
de los Reyes et al. 2015). 

3. K-band selection; The K-luminosity of galaxies mea¬ 
sures mass in old stars and therefore is a robust esti¬ 
mator of galaxy stellar mass (Broadhurst et al. 1992; 
Brinchmann & Ellis 2000). In this case, galaxies are se¬ 
lected by their luminosity in this band and then the red- 
shift is estimated using Spectral Energy Distribution 
(SED) htting (Drory & Alvarez 2008). However, the K- 
band magnitude limit of surveys can restrict the stellar 
mass distribution at the low mass end. This is one of the 
drawbacks of a magnitude-limited survey (Reddy et al. 
2012). Drory et al. (2005) showed that galaxies can 
also be selected by mass using I-band, yielding mass 
functions that are consistent with K-band selection. 

Besides the above selection methods that are typically used 
to construct parent samples of galaxies, observers adopt fur¬ 
ther selection criteria to create a sub-population within the 
initial sample with the desired properties and exclude other 
objects. Eor example, color cuts (e.g. BzK or U-V vs V-K) 
can be applied to the parent sample of galaxies (e.g. a K- 
selected parent sample), so that only star forming objects are 
included (Daddi et al. 2007; Guo et al. 2013). In addition, it 
is typical to separate galaxies between red-dead and blue- 
star forming galaxies. Elbaz et al. (2007) used the rest-frame 
colour-magnitude diagram (g-band centered at 4825 A) to 
separate between the two populations and excluded the red 
objects. This selection is used to reject possible contamina¬ 
tions from emission line objects whose emission is not due to 
star-formation and typically excludes passive or highly star 
forming systems with large contents of dust. In contrast to 
the above methods, stellar mass selection includes, besides 
Star Eorming Galaxies (SEGs), objects with high contents of 
dust and/or high stellar masses (Kajisawa et al. 2010; Karim 
et al. 2011; Bauer et al. 2011). 

The hrst question that arises is; Do samples of galaxies that 
were obtained using different selection criteria and wave¬ 
lengths produce SER-M* relations that provide similar con¬ 
strains for models and theory? A comparison between dif¬ 
ferent redshift results (z ~ 4 - 7) would suggest that the ob¬ 
served SER-M* relation could be signihcantly affected by 
the technique used to sample the galaxies (Katsianis et al. 
2015). 

In addition, different authors use different methods 
to obtain the intrinsic SERs and dust corrections of the 


observed galaxies. The most common techniques are; 

1) Conversion of IRh-UV luminosities to SERs; 

2) SED htting (e.g. Bruzual & Chariot 2003) to a range of 
wavelengths to obtain dust corrections and/or SERs. 

Both techniques have their advantages and shortcom¬ 
ings. However, the second question that arises is; Do the two 
above different methods produce SER-M* relations that 
are consistent with each other? Interestingly, Kajisawa et al. 
(2010) and Bauer et al. (2011) demonstrated that the above 
different techniques produce different SER-M* relations 
for the same sample of galaxies at high redshifts. When 
the authors used the conversion of IR and UV luminosities, 
instead of SED htting, to estimate dust corrections, they 
obtained higher values of SER at a hxed stellar mass. In 
addition, Utomo et al. (2014) claimed that the UV and 
IR luminosities overestimate SERs compared to the SED 
SERs by more than 1 dex for galaxies with Specihc SER 
(sSER s SER/M*), log(SSER) < -lOyr '. However, for 
the young highest star-forming galaxies in their sample the 
two methods to derive the SERs were found broadly con¬ 
sistent. Similar results were supported by Eumagalli et al. 
(2014). Boquien et al. (2014) argued that SERs obtained 
from modeling that takes into account only EUV and U 
bands are overestimated. Einally, Hayward et al. (2014) 
noted that the SERs obtained from IR luminosities (e.g. 
Noeske et al. 2007; Daddi et al. 2007) can be artihcially 
high. Despite the fact that SED htting techniques imply 
SER-M* relations that are inconsistent with results that 
rely solely on IRh-UV-SER conversions, it is quite common 
for compilation studies to combine SER(M*) that were 
obtained employing different methods (e.g. Behroozi et al. 
2013). In Table 1 we present a summary of the different 
observations used for this work. We include the technique, 
sample selection, mass limit and area of sky covered. We 
will make a brief comparison between these to investigate if 
there is a tension between the relations reported by different 
authors at z ~ 1 - 4. However, we will mostly focus on 
the SER-M* relations that simulations produce and how 
different they are from observations. 

Motivated by its importance in understanding galaxy evo¬ 
lution, a number of authors (e.g. Dave 2008; Dutton et al. 
2010; Einlator et al. 2011; Dayal & Eerrara 2012; Kannan 
et al. 2014; Sparre et al. 2015; Eurlong et al. 2015; Katsianis 
et al. 2015) have used hydrodynamic and semi-analytic mod¬ 
els to predict the SER-M* relation and its evolution. Numer¬ 
ical results (e.g. Dave 2008) predict a steeper relation than 
is found in observations. Eor example, the SER-M* rela¬ 
tion presented by Kannan et al. (2014) provides good agree¬ 
ment with the observed relation of Kajisawa et al. (2010) for 
z ~ 3. On the other hand, for redshift z ~ 2.2 the simulated 
galaxies have a star formation rate that is only half as large 
as observed. Katsianis et al. (2015) reported a discrepancy 
between the simulated SER-M* relation and the observed 
relations of Bouwens et al. (2012) and Heinis et al. (2014) 
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Table 1 Summary of the different observations used for this work. 


Publication 

Redshift 

Technique to obtain 
SFR-dust corrections 

Parent sample 
Wavelength/technique 

Further selection- 
Final sample 

Original IMF 

Noeske et al. (2007) 

0.8, 1.0 

EL-rIR 

Optical 

Main Sequence-SFG 

Kroupa (2001) 

Elbaz et al. (2007) 

1.0 

UV-rIR 

UV+optical 

Blue-SFG 

Kroupa (2001) 

Daddi et al. (2007) ‘ 

2.0 

UV-hIR 

K-hand 

BzK-SEG 

Salpeter(1955) 

Drory & Alvarez (2008) 

2.6-3.8 

SED(2800A)-As£d “ 

I-band 

None-Mag limited 

Salpeter (1955) 

Magdis et al. (2010) 

3.0 

UV-rIR 

Optical 

LBG 

Salpeter (1955) 

Kajisawa et al. (2010) 

0.75-3.0 

SED(2800A)-E(B-V)( p„-l-IR 

K-hand 

Mass 

Salpeter (1955) 

Kajisawa et al. (2010)^ 

0.75-3.0 

IR-i-UV 

K-hand 

Mass 

Salpeter (1955) 

Karim et al. (2011) 

0.8-3 

Radio 

Optical 

Mass 

Chahrier (2003) 

Bauer et al. (2011) 

1.75-2.75 

SED(2800 A)-UV slopes ed‘ 

Multi-wavelength 

Mass 

Salpeter (1955) 

Bauer etal. (2011)^ 

1.75-2.75 

IR+UV 

Multi-wavelength 

only 24 fim detected 

Salpeter (1955) 

Reddy etal. (2012) 

2.0 

UV-rIR 

UV-LBG 

LBG, Malmquist bias 

Salpeter (1955) 

Reddy etal. (2012)^ 

2.0 

UV-rIR 

UV-LBG 

LBG, bias corrected 

Salpeter (1955) 

Bouwens et al. (2012) 

3.8 

UV-UVslope 

UV-LBG 

LBG 

Salpeter (1955) 

Kashino et al. (2013) 

1.5 

Ha-l-UV 

Optical 

Bzk-SFG 

Salpeter (1955) 

Behroozi et al. (2013) 

0-4 

Comp 

Comp 

Comp 

Chahrier (2003) 

Bauer et al. (2013) 

Local 

Hcr-Balmer decrement 

Multi-wavelength 

Mass 

Chahrier (2003) 

Koyama et al. (2013) 

0.8, 2.2 

Ha-Garn & Best (2010) 

Ha 

Ha 

Salpeter (1955) 

Guo etal. (2013) 

0.7 

UV-rIR 

Multi-wavelength 

U-V vs V-k -SEG 

Chahrier (2003) 

Heinis et al. (2014) 

1.5, 3.0, 4.0 

UV-rIR 

UV-LBG 

LBG 

Chahrier (2003) 

Whitaker et al. (2014) 

0.5-2.5 

UV-hIR 

NIR 

U-V vs V-J -SEG 

Chahrier (2003) 

de Barros et al. (2014) 

3.0, 4.0 

SED-As£d " 

UV-LBG 

LBG 

Salpeter (1955) 

de los Reyes et al. (2015) 

0.8 

^a-AsED ^ 

Ha 

Ha 

Chahrier (2003) 

Salmon et al. (2015) 

4.0 

SED-As£d > 

Multi-wavelength 

None 

Salpeter (1955) 

Tomczak et al. (2016) 

0.8-4.0 

UV-rIR 

Multi-wavelength 

Mass 

Chahrier (2003) 


Notes: Column 1, publication reference of the observed SFR-M* relation; column 2, redshift used; column 3, technique 
and type of luminosity used to obtain the intrinsic SFR and dust corrections. These are: Emission lines (EL), Infrared (IR) 
luminosity, Ultra-Violet (UV) luminosity. Spectral Energy Distribution (SED) and radio luminosity. We include the Com¬ 
pilation (Comp) of studies presented in Behroozi et al. (2013); column 4, parent sample selection method/wavelength; 
column 5, final sample selection method. These are: Star Eormation Galaxy (SEG), mass, Ha and Lyman Break Galaxy 
(LBG) selection; column 6, mass limit of the observations; column 7, area covered by the survey in deg^', column 8, orig¬ 
inal IME adopted. We present the relations suggested by the authors corrected to a Chabrier (2003) IME when necessary 
(this conversion does not significantly affect the relations given since a similar change to SER and M* is used) along with 
results from cosmological hydrodynamic simulations in Eigure 2 (z ~ 3.1 and z ~ 3.8), Eigure 3 (z ~ 2.2 and z ~ 2.6), 
Eigure 4 (z ~ 2.0 and z ~ 1.5) and Eigure 5 (z ~ 0.8 - 1.15). 1) Daddi et al. (2007) used a fraction of the GOODS-N and 
GOODS-S fields, but did not report an exact area. 2) These observations are for the same sample of galaxies of the original 
work but adopt different methods for the determination of intrinsic SERs and dust corrections. 3) This set of observations 
is corrected for the effects of the Malmquist bias. 

a) SERs from rest-frame UV (2800 A) which was calculated from the galaxy SEDs (U, B, g, R, I, 834 nm, z, j and K 
bands). Dust corrections using the extinction curve of Calzetti (1997). 

b) SERs from rest-frame UV (2800 A) which was calculated from SEDs using multi-band photometry (U, B, V, i, z, J, H, 
K, 3.6 gtm, 4.5 //m, and 5.8 //m). 

c) SERs from Optical ACS/zgso-band. Dust corrections from UV slope calculated from multiwavelength SED-fitting. 

d) SERs from SEDs (B, V, I, z, U, R, J, H and K bands)H-the effect of nebular emission lines. Dust corrections using the 
extinction curve of Calzetti et al. (2000). 

e) SER from Ha luminosity. Dust corrections from the SEDs of galaxies at rest-frame UV and optical bands. 

Dust corrections using Pei (1992) and the extinction curve of 


the authors are broadly in agreement with the relation given 
by the compilation of observations used by Behroozi et al. 
(2013). However, the normalization of the simulated relation 
is significantly lower than the observational constraints at 
z ~ 1 and z ~ 2. In addition, the simulated relation is steeper 
than observed at all redshifts. In another study, Eurlong et al. 
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f) SFRs from SEDs (B, V, i, I, z, Y, J, JH and H bands). 
Calzetti et al. (2000). 

at z ~ 4. On the other hand, good agreement was found be¬ 
tween numerical results and the observations of Drory & 
Alvarez (2008) and de Barros et al. (2014). More recently, 
Sparre et al. (2015) used high resolution simulations to in¬ 
vestigate the SER-M* relation for redshifts z ~ 0 - 4 as part 
of the Illustris project (Genel et al. 2014). At z ~ 4 and z ~ 0, 
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(2015) showed agreement between simulated and observed 
specific SFR - stellar mass relations at z ~ 0 using high res¬ 
olution simulations from the EAGLE project (Schaye et al. 
2015). However, once again, the observed relations have a 
significantly higher normalization at z ~ 1 and z ~ 2. The 
tension between observations and simulations, especially at 
intermediate redshifts, implies either that the current models 
of galaxy evolution are incomplete or that observations are 
being misinterpreted. 

This paper is the third of a series, in which we seek to 
study a range of models for galaxy formation in hydrody- 
namical simulations through comparisons with observations 
of SER and stellar mass across cosmic time. In the first pa¬ 
per (Tescari et al. 2014) we constrained and compared our 
hydrodynamic simulations with observations of the cosmic 
star formation rate density and Star Eormation Rate Eunction 
(SERE) at z ~ 4 - 7. In the second paper (Katsianis et al. 
2015) we demonstrated that the same cosmological hydro- 
dynamic simulations reproduce the observed Galaxy Stellar 
Mass Eunction (GSME) and SER-M* relation for the same 
redshift interval. In this work, we extend the analysis pre¬ 
sented in Katsianis et al. (2015) down to z ~ 1, and criti¬ 
cally address the comparison of SER-M* relations with ob¬ 
servations obtained by using different analysis techniques. 
In section 2 we briefly summarize our numerical methodol¬ 
ogy. In section 3 we compare our results with the simulated 
SER-M* relations from different groups and demostrate the 
excellent consistency between numerical results from vari¬ 
ous projects (e.g. ANGUS, Illustris, EAGLE). In section 4 
we present the evolution of the SER-M* relation alongside 
observations from different groups that used different tech¬ 
niques to obtain their results. We draw our conclusions in 
section 5. In Appendix A we discuss how the uncertainty of 
the observed SER-M* relation has affected the comparison 
with simulations in the past. 

2 Simulations 

In this work we use the set of AustraliaN GADGET-3 early 
Universe Simulations (ANGUS) described in Tescari et al. 
(2014)'. We run these simulations using the hydrodynamic 
code P-GADGET3(XXL). We assume a flat A cold dark mat¬ 
ter (ACDM) model with Qom = 0.272, Qob = 0.0456, = 

0.728, n, = 0.963, Hq = 70.4 km s^' Mpc^' (i.e. h = 0.704) 
and cTg = 0.809. Our configurations have box size L — 2A 
Mpc/fi, initial mass of the gas particles Mgas = 7.32 x 10® 
Mo/fi and a total number of particles equal to 2 x 288^. All 
the simulations start at z = 60 and were stopped at z = 0 . 8 . 
Eor this work we use a Chabrier (2003) Initial Mass Eunction 
(IME) for all configurations. 

We explore different feedback prescriptions, in order to 
understand the origin of the difference between observed and 


*The features of our code are extensively described in Tescari et al. 
(2014) and Katsianis et al. (2015), therefore we refer the reader to those 
papers for additional information. 


simulated relationships. We also study the effects of metal 
cooling. We do not explore the broadest possible range of 
simulations, but concentrate on the simulations that can de¬ 
scribe the high-z star formation rate function and galaxy 
stellar mass function. We performed resolution tests for 
high redshifts (z ~ 4 - 7) in the appendix of Katsianis et al. 
(2015) and showed that our results converge for objects with 
logio < (M*/M 0 ) > 8.5. In Table 2 we summarise the main 
parameters of the cosmological simulations performed for 
this work. 


2.1 SNe feedback 


We investigate the effect of three different galactic winds 
schemes in the simulated SER-M* relation. We use the 
implementation of galactic winds of Springel & Hern- 
quist (2003). We assume the wind mass loading factor rj — 
M„/M* = 2 and a fixed wind velocity Vw = 450 km/s. In 
addition, we explore the effects of variable winds. We use 
a momentum driven wind model in which the velocity of 
the winds is proportional to the circular velocity Vcirc of the 
galaxy: 


Vw — 2 


GMhalo 


R 


— 2 X Vcirc. 


200 


and the loading factor rj. 


7] — 2x 


450 km/s 


( 1 ) 

( 2 ) 


where Mhaio is the halo mass and R 200 is the radius within 
which a density 200 times the mean density of the Universe 
at redshift z is enclosed (Barai et al. 2013). Eurthermore, 
we investigate the effect of the energy driven winds used by 
Puchwein & Springel (2013). In this case the loading factor 
is 


j] -2x 


1450 km/s 

\ Vw 


(3) 


while Vw = 2 X Vcirc 


2.2 AGN feedback 

In our scheme for Active Galactic Nuclei (AGN) feedback, 
when a dark matter halo reaches a mass above a given 
mass threshold M* = 2.9 x lO*" Mo/fi for the first time, it 
is seeded with a central Super-Massive Black Hole (SMBH) 
of mass Mseed = 5.8 X 10^ Mo/fi (provided it contains a min¬ 
imum mass fraction in stars /* = 2.0 x 10“"'). Each SMBH 
will then increase its mass by accreting local gas from a max¬ 
imum accretion radius R^c - 200 kpc/fi. In this scheme we 
allow the presence of a black hole in low mass halos, and 
at early times. The AGN feedback prescription that we use 
combined with efficient winds is successful at reproducing 
the observed SERE (Tescari et al. 2014) and GSME (Katsia¬ 
nis et al. 2015) for redshifts 4 < z < 7. 
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Table 2 Summary of the different runs used in this work. 


Run 

IMF 

Box Size 
[Mpc/Zi] 

Ntot 

Mdm 

[Mo//i] 

Mgas 

[Mo/M 

Comoving Softening 
[kpc/h] 

Feedback 

Ch24.eA_CsW 

Chabrier 

24 

2 X 288-^ 

3.64x10' 

7.32 X lO" 

4.0 

Early AGN + Constant strong Winds 

Ch24^AjjW 

Chabrier 

24 

2 X 288-^ 

3.64x10' 

7.32 X lO" 

4.0 

Early AGN + no Winds 

Ch24JXF 

Chabrier 

24 

2 X 288" 

3.64x10' 

7.32 X 10” 

4.0 

No Feedback 

Ch24.eAMDW‘‘ 

Chabrier 

24 

2 X 288-' 

3.64x10' 

7.32 X 10*' 

4.0 

Early AGN + 
Momentum-Driven Winds 

Ch24j;Aj:DW‘’ 

Chabrier 

24 

2 X 288-^ 

3.64x10' 

7.32 X lO" 

4.0 

Early AGN -t- 
Energy-Driven Winds 

Ch24JZ:c^A.CsW‘ 

Chabrier 

24 

2 X 288" 

3.64x10' 

7.32 X 10” 

4.0 

Early AGN -i- Constant strong Winds 
Metal cooling 

Ch24JZc_eAj:DW^ 

Chabrier 

24 

2 X 288" 

3.64x10' 

7.32 X lO” 

4.0 

Early AGN -i- Energy-Driven Winds 
Metal cooling 


Notes: Column 1, run name; column 2, Initial Mass Function (IMF) chosen; column 3, box size in comoving Mpc//i; 
column 4, total number of particles (Ntot = Ngas + Ndm); column 5, mass of the dark matter particles; column 6, initial 
mass of the gas particles; column 7, Plummer-equivalent comoving gravitational softening length; column 8, type of 
feedback implemented. See section 2 and Tescari et al. (2014) for more details on the parameters used for the different 
feedback recipes, (a): in this simulation we adopt variable momentum-driven galactic winds (Subection 2.1). {b): in this 
simulation we adopt variable energy-driven galactic winds (Subsection 2.1). (c): in these simulations the effect of metal¬ 
line cooling is included (Subsection 2.3). For all the other runs we use cooling tables for gas of primordial composition 
(H H- He). 


2.3 Metal cooling 

Our code follows the evolution of 11 elements (H, He, C, 
Ca, O, Ne, Mg, S, Si and Fe) released from supernovae 
(SNIa and SNII) and low and intermediate mass stars self- 
consistently (Tornatore et al. 2007). Radiative heating and 
cooling processes are included following Wiersma et al. 
(2009). We assume a mean background radiation composed 
of the cosmic microwave background and the Haardt & 
Madau (2001) ultraviolet/X-ray background from quasars 
and galaxies. Contributions to cooling from each one of the 
eleven elements mentioned above have been pre-computed 
using the Cloudy photo-ionisation code (last described in Fer- 
land et al. 2013) for an optically thin gas in (photo)ionisation 
equilibrium. In this work we use cooling tables for gas of 
primordial composition (H + He) as the reference configura¬ 
tion. To test the effect of metal-line cooling, we include it in 
two simulations (Ch24JZcueA_EDW and Ch24JZc-eAJTsW). 


3 Comparison between SFR-M* relations from 
different simulations 

Cosmological hydrodynamic simulations provide a power¬ 
ful tool to investigate and predict properties of galaxies and 
their distribution. Recently, the Illustris (Sparre et al. 2015) 
and EAGLE (Eurlong et al. 2015) projects have used high 
resolution simulations and tried to reproduce the observed 
SER-M*. In addition, semi-analytic models have been 
used to reproduce the observed relations at various redshifts 
(Dutton et al. 2010). The main questions that arise are: 


a) Are simulations and theory capable of reproducing the 
observed SER-M* relations? 

b) Are the results of simulations that are tuned to reproduce 
certain observables consistent with each other? 

c) Cosmological hydrodynamic simulations have been 
evolving remarkably in the last decade in terms of resolution 
and box size. However, numerical modeling of the inter¬ 
stellar medium (ISM) and star formation physics remained 
essentially the same (e.g. stochastic formation of star 
particles in a multiphase ISM). Therefore, how different are 
the SER-M* relations found by state-of-the-art simulations 
with respect to those found in the past? 

In Eigure 1 we present a comparison between the sim¬ 
ulated median SER-M* relations of our reference model 
(Project name: ANGUS, black solid line), and those pre¬ 
sented in Dutton et al. (2010, blue dotted line, semi-analytic 
model), Sparre et al. (2015, Illustris, magenta dashed line) 
and Eurlong et al. (2015, EAGLE, red dashed line). Eor this 
analysis we do not include galaxies that have masses lower 
than the confidence limit of 10® M© from our simulations. 
This is done to make a meaningful comparison with the Il¬ 
lustris and EAGLE projects. We see the agreement between 
different groups is excellent despite the fact that they used 
different resolutions and box sizes (in Table 3 we include 
some details for the runs that were used to produce the sim¬ 
ulated SER-M* shown in Eigure 1). Models that are tuned 
to reproduce certain observables (eg. ANGUS - cosmic star 
formation rate density evolution, EAGLE - GSME at z ~ 0, 
Dutton et al. (2010) - SER-M* relation at z ~ 0) produce 
similar results for the star formation rate main sequence. It is 
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Figure 1. Median values of the SFR-M* relations from different cosmological hydrodynamic simulations for z ~ 0 - 4. The black line is the median line 
through all the points of the scatter plot for our reference model {Ch24.eAJLDW). The blue dotted line is the median fit of the scatter plot obtained with 
the Semi-Analytic Model (SAM) of Dutton et al. (2010). The magenta dashed line is the median line of the scatter plot presented by Sparre et al. (2015, 
Illustris project). The red dashed line is the median line of the scatter plot presented hy Furlong et al. (2015, EAGLE project). We cut our SFR(M*) under 
our confidence mass limit of 10® Mq to make a meaningful comparison with the Illustris and EAGLE projects. There is an excellent agreement between the 
results from cosmological hydrodynamic simulations run by different groups. At each redshift, a panel showing ratios between the different simulations and 
observations with the Ch24.eAJIDW (black solid line) is included. 


worth to mention that even under our mass confidence limit 
our results are in agreement with Furlong et al. (2015). 

Dave (2008) using GADGET-2 (Springel 2005) demon¬ 
strated that simulations of galaxy formation produce sim¬ 
ilar SFR-M* relations, to a large extent independently of 
modeling details (e.g feedback prescriptions). According to 
the authors, this is a generic consequence of smooth and 
steady cold accretion which dominates the growth of the 
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simulated objects. In addition, Dutton et al. (2010) state that 
the SFR-M* relation is generally found to be independent 
of feedback, since feedback regulates the outflow rate, and 
mostly acts to shift galaxies along the SFR sequence, leav¬ 
ing the zero point of the relation invariant. We display the ex¬ 
cellent agreement between the SFR-M* relations found in 
different cosmological hydrodynamic simulations (ANGUS, 
Illustris, EAGLE) and the semi-analytic results of Dutton 
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Table 3 Summary of the different simulated SFR—M* relations used for Figure 1. 


Publication 

redshift 

Box Size 
[Mpc/h] 

Ntot 

Mdm 

[Mo/A] 

Mgas 

[Mo/A] 

Comoving Softening 
[kpc//i] 

Feedback 

Ch24jeAMDW (this work) 

0, 1.15,2.0,3.8 

24 

2 X 288-’ 

3.64x10' 

7.32 X 10” 

4.0 

Early AGN + 
Energy-Driven Winds 

Sparre et al. (2015) 

0, 1,2,4 

75 

2 X 1820’ 

4.41x10” 

8.87x10” 

0.704 

AGN + Stellar 
Vogelsberger et al. (2013) 

Furlong et al. (2015) 

0, 1,2 

70.4 

2 X 1504’ 

6.83x10” 

1.27x10” 

1.87 

AGN + Stellar 

Crain etal. (2015) 

Dutton et al. (2010) 
Semi-analytic 

0, 1,2,4 

NA 

NA 

NA 

NA 

NA 

Stellar 

Dutton et al. (2010) 


Notes: Column 1, publication reference; column 2, redshifts considered; column 3, box size of the simulation in comoving 
Mpc/h; column 4, number of particles used; column 5, mass of the dark matter particles; column 6, initial mass of the gas 
particles; column 7, comoving softening length; column 8, feedback prescriptions used. The box size and masses of the 
dark matter and gas particles are in Mpc//z and respectively, rescaled to our adopted cosmology {h = 0.704). 


et al. (2010). This strongly suggests that the slope and nor¬ 
malization of the relation have their origins in fundamental 
principles and assumptions commonly adopted in numerical 
codes, while the small differences in feedback prescriptions 
play a negligible role. Our results for the local Universe are 
consistent with the relations of other groups, despite the fact 
that our box size and resolution are not sufficient to robustly 
probe other properties of galaxies at z ~ 0. 

Increasing the resolution provides a description of smaller 
masses and scales. In addition, decreasing the softening 
length will better resolve the feedback mechanisms. For this 
reason, some authors claim that an exact convergence when 
resolution is changed is not to be expected (Schaye et al. 
2015). However, Murante et al. (2015), who used simu¬ 
lations of disc galaxies (based on the GADGET-3 code), 
demonstrated that numerical results are remarkably stable 
against resolution even in runs dedicated to galactic scales. 
For their ISM multiphase model the authors employed dif¬ 
ferential equations that describe the evolution of a system 
composed of cold clouds (where stars form) embedded in 
hot ambient gas, at unresolved scales (MUPPI). The results 
from different runs are stable as resolution is decreased even 
by a factor of 8. In particular, morphology-related quanti¬ 
ties, such as rotation curves and circularity histograms, vary 
by less than 10 per cent. The SFR varies approximately by 
the same amount. Murante et al. (2015) suggest that reduc¬ 
ing the softening in simulations of disk galaxies by a factor 
of 6 induces effects related to numerical heating that change 
their morphologies and central velocities. Doubling the soft¬ 
ening parameter results in thicker and less extended discs 
and increases the bulge mass. However, the results for the 
integrated properties are not significantly affected. The re¬ 
markable agreement between the simulated SFR - M* rela¬ 
tions presented in Figure 1 suggests that cosmological sim¬ 
ulations are able to produce realistic populations of galax¬ 
ies within representative cosmological volumes, even at rela¬ 
tively modest resolutions. Moreover, it is a strong indication 
that numerical codes that rely on similar multiphase models 


and are tuned to reproduce statistical properties of galaxies 
are “bound” to produce similar results by construction, re¬ 
gardless of resolution and box size. 

In the following sections we investigate in more detail the 
redshift evolution of the simulated SFR-M* relations in the 
Angus project and critically compare the numerical results 
with observations. 


4 Evolution of the SFR-M* relation from z ~ 4 to 

z ~ 1 

In Katsianis et al. (2015) we demonstrated that different ob¬ 
servations of the SFR-M* relation are in tension for z ~ 
4-7. This discrepancy was attributed to the different se¬ 
lection methods and techniques for the determination of the 
intrinsic SFRs and dust corrections of the observed galax¬ 
ies. In this section we extend that work to z ~ 1 - 4, in or¬ 
der to investigate if observations in this redshift interval are 
also in tension, and address in more detail how sample se¬ 
lection and methodology can affect the comparison with cos¬ 
mological simulations. We present the SFR-M* relation ob¬ 
tained from our hydrodynamic simulations along with ob¬ 
servations from different groups in Figs. 2 (z ~ 3.1 - 3.8), 
3 (z ~ 2.2 - 2.6), 4 (z ~ 1.5 - 2.0) and 5 (z ~ 0.8 - 1.15). In 
the left panels of each figure we show the scatter plots of the 
SFR-M* relation for our reference model Ch24_eA_EDW 
(grey points). The black line is the median line through 
all points of the scatter plot. Our reference model, which 
combines a Chabrier IMF, early AGN feedback and vari¬ 
able energy driven winds, is able to reproduce the observed 
galaxy stellar mass function, star formation rate function 
and cosmic star formation rate density for z ~ 1 - 7 galax¬ 
ies (Tescari et al. 2014; Katsianis et al. 2015). In the right 
panels we compare the median lines of the scatter plots for 
all the runs presented in Table 2. The orange vertical line at 
10® Mo is the confidence limit of our simulations and the 
mass limit of most observations for the redshifts considered 
in this work. 
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Figure 2. Left panels: scatter plot of the SFR-M* relation for our fiducial run Ch24.eA^DW (grey points) at redshifts ^ ~ 3.1 and 3.8. In each panel, the 
black solid line is the median line through the points of the scatter plot. Overplotted are the observed galaxy SFR(M*) relations from Drory & Alvarez (2008, 
I-band selected sample, SFRs(SED) - brown dashed line), Magdis et al. (2010, Lyman-break selected sample, SFRs(UV+IR) - black dotted line), Kajisawa 
et al. (2010, mass-selected sample - the dark green stars represent SFRs that were obtained using UV+IR luminosities, while the blue stars were obtained 
using the SED fitting technique), Bouwens et al. (2012, Lyman-break selected sample, SFRs(UV+IRX-yS) - black triple dot-dashed line), Heinis et al. (2014, 
Lyman-break selected sample, SFRs(UV+IR) - orange triple dot-dashed line), de Barros et al. (2014, Lyman-break selected sample, SFRs(SED) - blue filled 
circles with error bars), Salmon et al. (2015, multi-wavelength derived redshifts, SFRs(SED) - reverse green triangles with error bars) and Tomczak et al. 
(2016, mass-selected sample, SFRs(UV+IR) - black open circles with error bars). Right panels: median lines of the SFR—M* scatter plots for all the runs of 
Table 2. In these panels we do not present the analytic expressions of the observed relations for the sake of clarity. 


The scatter plot of the simulated SFR-M* for our refer¬ 
ence model at z ~ 3.8 (top left panel of Figure 2) is consis¬ 
tent with the results of the I-band selected sample of Drory & 
Alvarez (2008), the Lyman-Break selected sample of de Bar¬ 
ros et al. (2014) and the multi-wave length results of Salmon 
et al. (2015). On the contrary, the SFR-M* relations ob¬ 
tained by Bouwens et al. (2012) and Heinis et al. (2014), 
who used the Lyman-Break technique, imply ~ 3-5 times 
higher SFRs at a fixed stellar mass. All the above authors 
use various methods to obtain the intrinsic SFRs and dust 


attenuation effects. Drory & Alvarez (2008) and de Barros 
et al. (2014) used SED fitting techniques. On the other hand, 
Bouwens et al. (2012) and Heinis et al. (2014) estimated the 
dust attenuation effects and SFRs in their sample using the 
IRX-yS relation (Meurer et al. 1999), stacking techniques and 
the Kennicutt (1998) relation. The tension between the above 
results may be due to the fact that the authors used different 
methods to obtain the intrinsic properties of galaxies and/or 
due to selection methods. The SFR-M* given by Tomczak 
et al. (2016) represents a multi-wavelength sample of galax- 
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ies, where the SFRs were obtained from converting UV+IR 
luminosities. The relation implies lower values of SFR at a 
fixed stellar mass than the Lyman break selected sample of 
Heinis et al. (2014) who used UV+IR-SFR conversions as 
well. Heinis et al. (2014) discussed the possible impact of 
the Lyman-break selection on their retrieved SFR-M* rela¬ 
tion, and state that their sample may not include a significant 
number of objects since they preferably select star forming 
objects. They point out that their shallow slope, ~ 0.7, might 
be also caused by the fact that they are selecting galaxies 
by their UV flux, which does not pick objects with low star 
formation rates at fixed stellar mass or dusty massive ob¬ 
jects. The above would suggest that the differences between 
the results of the two authors could be attributed to the dif¬ 
ferences of their selection methods. On the other hand, the 
results of Tomczak et al. (2016) suggest a higher normaliza¬ 
tion than the observations of Salmon et al. (2015). Both au¬ 
thors used multi-wavelength samples of galaxies but Salmon 
et al. (2015) used an SED method to obtain SFRs instead of 
a UVh- 1R conversion. As discussed in the introduction it is 
quite possible that different techniques produce different re¬ 
sult. Differences between Salmon et al. (2015) and Tomczak 
et al. (2016) maybe can be attributed to the technique used 
to obtain dust corrections and SFRs. For redshiftz = 3.8, our 
configurations that combine different feedback schemes re¬ 
semble each other closely (top right panel of Figure 2). This 
follows the results of Katsianis et al. (2015), who showed 
that different feedback prescriptions result in roughly the 
same SFR-M* relation for z ~ 4 - 7. 

At z ~ 3.1, our reference model (Ch24_eA_EDW, bottom 
left panel of Figure 2) is consistent with the the mass se¬ 
lected sample of Kajisawa et al. (2010), and the Lyman- 
break selected sample of de Barros et al. (2014). Kajisawa 
et al. ( 2010 ) used two different methods to obtain the in¬ 
trinsic SFRs at a fixed stellar mass for the same sample of 
galaxies. In the first case, they used the sum of IR and UV 
light to estimate the observed SFRs for objects that were de¬ 
tected having an IR 24 fim flux (originating from dusty high 
star forming galaxies). For the other objects (undetected at 
24 yum) the authors used SED fitting techniques and UV lu¬ 
minosities. The dark green stars of Eigure 2 are obtained 
with this methodology. In the second case, SED fitting tech¬ 
nique and UV luminosities were used for the whole sample 
of galaxies. They find good agreement between both meth¬ 
ods at z ~ 3.0, even though the SERs from the combination 
of UV and IR light are higher than those found by SED fit¬ 
ting. Kajisawa et al. (2010) stressed that, at z ~ 3, the ratio 
log(SER/s+[/v/SER[/v^go) for galaxies with 24 fim flux is as 
high as ~ 0.63. This is an example of how methodology can 
affect the determination of SERs in a sample of galaxies. The 
relations given by the Lyman-break selection of Magdis et al. 
(2010) and Heinis et al. (2014) have significantly higher nor¬ 
malizations than those found by SED measurements, mass 
selected samples and cosmological simulations. 

We see the effect of metal cooling and different feed¬ 
back prescriptions among the different simulations at z ~ 


3.1 in the bottom right panel of Eigure 2. The Ch24J>}F 
(no feedback) and Ch24jeAjiW (early AGN, no winds) 
runs are almost identical. This means that the effect of our 
AGN feedback model on the simulated SER-M* relation 
is small. Eurthermore, we can compare the Ch24_eAjiW 
and Ch24jeAjCsW (early AGN, constant strong winds) to 
gain insight into how the constant energy wind model of 
Springel & Hernquist (2003) affects the simulated relation. 
We see that the star formation rate at fixed mass for ob¬ 
jects with stellar mass logjQ(M*/Mo) < 10.0 is lower for 
the Ch24jeAjCsW run. Above the mass limit of 10® M©, 
the Ch24_eAJADW (early AGN, momentum-driven winds) 
and Ch24^A_EDW models are consistent with the config¬ 
urations that have no winds. This means that the effect of 
our variable wind models do not change the slope of the 
SER-M* relation. The only small difference between our 
runs is found for objects with log[Q(M*/M 0 ) < 8.5, where 
there are no observations to constrain the results and feed¬ 
back is not well resolved. Note that Dave (2008) and Dut¬ 
ton et al. ( 2010 ) also suggest that runs with different feed¬ 
back prescriptions results in similar SER-M* relations^. Eor 
this work we use a set of physically plausible cases that can 
produce realistic star formation rate and stellar mass func¬ 
tions in our simulations (Tescari et al. 2014; Katsianis et al. 
2015). In addition, by comparing the Ch24JZc_eAjCsW and 
Ch24_eAJ2sW we see that metal cooling does not signifi¬ 
cantly change the simulated SER-M* relation. This is due 
to the fact that when metal cooling is included gas can cool 
more efficiently. As a result the SER increases and, corre¬ 
spondingly, the stellar mass increases moving galaxies along 
the SER-M* relation without affecting it considerably. 

The scatter plot of the simulated SER-M* relation for 
our fiducial model at z ~ 2.6 (top left panel of Eigure 3) is 
consistent with the results of Drory & Alvarez (2008) and 
Bauer et al. (2011). The open magenta triangles show the 
median SER that relied on adding IR and UV luminosities 
(SER/s+(/v) for the sample of galaxies of Bauer et al. ( 2011 ), 
which were detected only at 24 fim. The authors state that, at 
z > 2.5, the SER/«+(/v is greater than the SER obtained from 
SED and UV light by an average factor of 5. The full mass 
selected sample of Bauer et al. (2011) with dust correction 
laws that rely on SED fitting is almost 0.7 dex lower and 
in excellent agreement with numerical results (filled red tri¬ 
angles). This comparison points out how selection and dust 
correction methods can affect the relation reported by ob¬ 
servers at z ~ 2.6. Bauer et al. (2011) noted that the best way 
to robustly determine the SER and the amount of dust extinc¬ 
tion for each galaxy is to calculate the ultraviolet slope SED- 
fitting. We find that cosmological hydrodynamic simulations 
from various groups have a good consistency with these SED 
measurements of the SER-M* relation. On the contrary, 
simulations are unable to reproduce the results of Karim 


^Extreme feedback recipes can shape the SFR-M* relation (Haas et al. 
2013). However, these runs are unable to produce galaxies with realistic 
SFRs and stellar masses. 
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Figure 3. Left panels: scatter plot of the SFR-M* relation for our fiducial run Ch24.eA^DW (grey points) at redshifts ^ ~ 2.2 and 2.6. In each panel, the 
black solid line is the median line through the points of the scatter plot. Overplotted are the observed galaxy SFR(M*) relations from Drory & Alvarez (2008, 
I-band selected sample, SFRs(SED) - brown dashed line), Karim et al. (2011, SFG/mass-selected sample, SFRs(radio) - magenta/black circles with error 
bars), Bauer et al. (2011, mass-selected sample/SFRs that are obtained from SED fitting - red triangles), Bauer et al. (2011, 24 fim selected sample/SFRs 
that are obtained from UV+IR luminosities - magenta open triangles), Reddy et al. (2012, Lyman-break selected sample, SFRs(UV+IR) - dark green open 
squares are not corrected for incompleteness, dark green filled squares are the con'ected results), Koyama et al. (2013, Hn-selected sample, SFRs(HQf) - 
orange solid line), Whitaker et al. (2014, U-V vs V-J -SFG sample, SFRs(UV+IR) - orange open diamonds) and Tomczak et al. (2016, mass-selected sample, 
SFRs(UV+IR) - black open circles with error bars). Right panels: median lines of the SFR-M* scatter plots for all the runs of Table 2. In these panels we 
do not present the analytic expressions of the observed relations for the sake of clarity. 


et al. (2011), shown by the magenta filled circles (SFGs) 
and black filled circles (mass selected sample), who used 
radio luminosities to obtain the intrinsic SFRs. Speagle et al. 
(2014) noted that the SFRs obtained from radio luminosities 
are overestimated and in tension with other SFR indicators. 
We see that the SED observations and numerical results have 
SFRs at a fixed stellar mass lower by a factor of ~ 6. For red- 
shift z = 2.6, the Ch24jeA_CsW and Ch24JZc^A_CsW runs 
underpredict the SFR at a fixed M* for objects with stellar 
masses logjQ(M*/Mo) < 10.2. The rest of the configurations 


are consistent with each other and with the mass selected 
observations of Bauer et al. (2011, top right panel of Figure 
3). By comparing the Ch24JZc_eAJi,DW and Ch24^AJ2DW 
we see that metal cooling does not affect the simulated 
SFR-M* relation when energy driven winds are used. This 
is true for all redshifts considered in this work. 

The simulated SFR-M* relation for z ~ 2.2 (bottom left 
panel of Figure 3) is consistent with the I-band selected sam¬ 
ple of Drory & Alvarez (2008) and mass selected sample of 
Bauer et al. (2011). The results of Bauer et al. (2011) that 
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were obtained using two different dust correction methods 
(open magenta triangles and red filled triangels) agree bet¬ 
ter at z < 2.25. In the bottom left panel of figure 3 we also 
present the LEG observations of Reddy et al. (2012). We 
emphasize the comparison between the biased and unbiased 
results of Reddy et al. (2012) who investigated a set of galax¬ 
ies at redshifts 1.5 < z < 2.6. The authors quantified the ef¬ 
fects of the non detection of faint objects in their flux lim¬ 
ited selection. We see that the effect of the Malmquist bias 
(preferential selection of the most luminous-SFR galaxies at 
a fixed stellar mass) is important for low mass galaxies, and 
that the slope of the biased SFR-M* relation is therefore 
artificially shallow. The correction makes the slope steeper 
with an exponent close to unity, something that is in accor¬ 
dance with the predictions from cosmological simulations. 
We note that the constrains that include the correction for 
the Malmquist bias are in very good agreement with the up¬ 
dated observations of Whitaker et al. (2014) and Tomczak 
et al. (2016), which are able to better probe low mass objects. 
Reddy et al. (2012), Whitaker et al. (2014) and Tomczak 
et al. (2016) used combinations of IR and UV luminosities to 
obtain the intrinsic SFRs and dust corrections at z ~ 2.2 and 
maybe this is the reason why they predict larger normaliza¬ 
tions for their SFR-M* relation with respect to SED obser¬ 
vations and simulations. The Ha selected sample of Koyama 
et al. (2013) implies higher values of SFR at a fixed stellar 
mass than the mass selected sample of Bauer et al. (2011), 
and this maybe is due to the fact that it preferably samples 
high star forming galaxies ( Ha selection). Radio SFRs are 
almost a factor of ~ 4 larger than simulations estimates at 
z ~ 2.2.Forredshiftz = 2.2, the configurations with constant 
energy driven winds underpredict the SFR at a fixed mass 
for objects with log[o(M*/M0) < 10.3 (bottom right panel 
of Figure 3). However, the simulation with metal cooling and 
constant energy driven winds slightly overpredicts the SFR 
at fixed stellar mass for objects with log[Q(M*/M0) > 10.5 . 

The SFR-M* relation for our reference model at z ~ 2.0 
(top left panel of Figure 4) is consistent with the mass se¬ 
lected samples of Kajisawa et al. (2010, SED) and Bauer 
et al. (2011). However, we note that the observed relations 
are shallower. In contrast, the results from the BzK-SEGs of 
Daddi et al. (2007), the mass selected sample of Kajisawa 
et al. (2010, UVh-IR) and the U-V vs V-J -SEGs of Whitaker 
et al. (2014) imply a significantly higher normalization than 
SED measurements. Bauer et al. (2011) state that they find a 
flattened relation relative to Daddi et al. (2007). According to 
Bauer et al. (2011) this is either due to the fact that they are 
using a mass-complete sample instead of just star-forming 
galaxies, or the overestimation of the dust correction applied 
by Daddi et al. (2007) who used a combination of IR and UV 
luminosities instead of a SED analysis. Eurthermore, Hay¬ 
ward et al. (2014) note that the overestimation of SERs from 
IR luminosities may have played an important role in the ob¬ 
served SER-M* relations of Daddi et al. (2007) and suggest 
that the methodology used by the authors may have overes¬ 
timated the SERs at a fixed stellar mass. We see that the SER 


selected results of Daddi et al. (2007, red triple dot-dashed 
line - IRh-UV) and the mass selected sample of Kajisawa 
et al. (2010, dark green symbols - UVh-IR) are in excellent 
agreement. This points to the direction that the BzK-SEG 
selection of Daddi et al. (2007) did not considerably affect 
the SER-M* relation and the methodology for dust correc¬ 
tions and determinations of intrinsic SERs is mostly respon¬ 
sible for the tension with the results of Kajisawa et al. (2010, 
blue symbols-UVs£o) and Bauer et al. (2011). Kajisawa 
et al. (2010) state that at z ~ 2 the ratio of SERuv+ir and 
SERuvsed is log(SERiR+uv/SERuvsED) = 0-37, for the same 
sample of galaxies. We also see from Eigure 4 that the dif¬ 
ference between the two methods for the determination of 
the intrinsic SERs at a fixed stellar mass is ~ 0.75 dex for 
objects with > 10^® ^ M©. Simulations are more consistent 
with SED observations at z ~ 2, while the use of IR light 
predicts larger SERs from both. Numerical results are con¬ 
sistent with the observations of Tomczak et al. (2016) for 
high mass objects > lO'® ® M©, but at lower masses the ob¬ 
servations suggest significantly higher SERs at a fixed stel¬ 
lar mass. Overall the observed relation is found to be shal¬ 
lower. Eor redshift z = 2.0, the simulations with constant en¬ 
ergy driven winds underpredict the SER at a fixed mass for 
objects with log[o(M*/M©) < 10 (top right panel of Eigure 
4). The simulation with metal cooling and constant energy 
driven winds overpredicts the SER at fixed stellar mass for 
objects with log[Q(M*/M0) > 10. 

The simulated SER-M* relation for the Ch24_eA_EDW at 
z ~ 1.5 (bottom left panel of Eigure 4) is steeper with lower 
normalization than the results of Heinis et al. (2014). On the 
other hand, the normalization of the SER-M* relation ob¬ 
tained by their sample could be possibly larger due to their 
Lyman-break selection and/or the methodology used to re¬ 
cover the intrinsic properties of galaxies. We see that the 
results from Lyman-break and SEGs are in excellent agree¬ 
ment (Kashino et al. 2013; Heinis et al. 2014; Whitaker et al. 
2014) and this points to the direction that they select similar 
high star forming systems. The agreement of the simulations 
with the radio SERs of Karim et al. (2011) is improved (es¬ 
pecially for the mass selected sample). 

The scatter plot of the simulated SER-M* relation for 
our reference model at z ~ 115 (top left panel of Eigure 5) 
is consistent with the results of the Kajisawa et al. (2010, 
for M* < 1O'*''®M0) and Karim et al. (2011). However, all 
the observations suggest significantly higher SERs at a fixed 
stellar mass for small objects. We see though that the ten¬ 
sion between the observed relations is much less. The dif¬ 
ference between Kajisawa et al. (2010, dark green symbols- 
UVh-IR) and Kajisawa et al. (2010, blue symbols-UVsBo) 
is much smaller than that found at higher redshifts. The 
authors state that log(SER/R+(/v/SERj/yj^^) is ~ 0.25 and 
~ 0.19 for z ~ 1.25 and z ~ 0.75, respectively. The differ¬ 
ence between the two methods to recover the intrinsic SER 
at a fixed stellar mass for objects with M* > 10'°^ M© is 
~ 0.35 dex and ~ 0.25 dex for z ~ 1.25 and z ~ 0.75, re¬ 
spectively (considerably lower than the 0.75 dex found at 
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Figure 4. Left panels: scatter plot of the SFR-M* relation for our fiducial run Ch24.eAJiDW (grey points) at redshifts z ~ 1.5 and 2.0. In each panel, 
the black solid line is the median line through the scatter plot. Overplotted are the observed galaxy SFR(M*) relations from Daddi et al. (2007, BzK- 
SFGs, SFRs(UV+IR) - red triple dot-dashed line), Kajisawa et al. (2010, mass-selected sample, SED/UV-hIR - blue stars/green stars), Karim et al. (2011, 
SFGs/mass-selected sample - magenta/black circles), Bauer et al. (2011, mass-selected sample, SFRs(SED) - red triangles), Kashino et al. (2013, Lyman- 
break selected sample, SFRs(UV-l-Hci') - black dashed line), Fleinis et al. (2014, Lyman-break selected sample, SFRs(UV-l-IR) - orange triple dot-dashed line), 
Whitaker et al. (2014, U-V vs V-J -SFG sample, SFRs(LrV-l-IR) - orange open diamonds) and Tomczak et al. (2016, mass-selected sample, SFRs(UV-hIR) - 
hlack open circles with en'or bars). Right panels: median lines of the star SFR—M* scatter plots for all the runs of Table 2. In these panels we do not present 
the analytic expressions of the observed SFR—M* relations for the sake of clarity. 


z ~ 2.0). The turn over of the relation found in the mass 
selected sample of Kajisawa et al. (2010) and the I-band 
selected sample of Drory & Alvarez (2008) could be at¬ 
tributed to the fact that the authors include massive quiescent 
galaxies (Drory & Alvarez 2008). However, recent observa¬ 
tions of SFGs that do not include quiescent objects show as 
well a curvature (Whitaker et al. 2014). Simulations (includ¬ 
ing Angus, Illustris and EAGLE) are unable to reproduce 
this behaviour. It is quite possible that current mechanisms 
implemented in state-of-the-art cosmological simulations to 
decrease the SERs of high mass objects (like e.g. AGN feed- 
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back) are not efficient enough. Once again, we see that the 
run with constant energy driven winds (Ch24jeA_CsW) un¬ 
derpredicts the SER at a fixed mass (top right panel of Figure 
5). The simulation with constant energy driven winds and 
metal cooling underpredicts the SFR at a fixed stellar mass 
for logio(M*/Mo) < 10.4. On the contrary, for objects with 
logio(M*/Mo) > 10.5 the SFR is overpredicted. The other 
configurations are consistent with each other. 

The simulated SFR-M* relation for our fiducial model at 
z ~ 0.8 (bottom left panel of Figure 5) is consistent with the 
mass selected samples of Kajisawa et al. (2010) and Karim 
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Figure 5. Left panels: scatter plot of the SFR-M* relation for our fiducial run Ch24.eA^DW (grey points) at redshifts z ~ 0.8 and 1.15. In each panel, the 
black solid line is the median line through the scatter plot. Overplotted are the observed galaxy SFR(M*) relations from Kajisawa et al. (2010, mass-selected 
sample - SED/blue stars, IR+UV/dark green stars), Karim et al. (2011, SFGs/mass-selected sample - magenta/black circles), Guo et al. (2013, U-V vs V-K 
SFGs, SFRs(UV+IR) - red dashed line), Koyama et al. (2013, Hcr-selected sample, SFRs(Ho') - orange solid line), Whitaker et al. (2014, U-V vs V-J -SFG 
sample, SFRs(UV+IR) - orange open diamonds) and Tomczak et al. (2016, mass-selected sample, SFRs(UV+IR) - black open circles with error bars). Right 
panels: median lines of the SFR-M* scatter plots for all the runs of Table 2. In these panels we do not present the analytic expressions of the obseiwed 
SFR-M* relations for the sake of claiity. 


et al. (2011), the Ha selected sample of Koyama et al. (2013) 
and the U-V - V-k selected sample of Guo et al. (2013). The 
updated relation given by Guo et al. (2013) is steeper than 
other estimates in the literature (e.g. Noeske et al. 2007). 
Guo et al. (2013) noted that taking into account low mass 
objects is critical for determining the slope of the SFR-M* 
relation and the reason why they find a steeper slope (with 
an exponential slope close to unity) is that they can take into 
account objects that previous surveys did not. If the authors 
fit only their relation to the 24 jum detected galaxies, they 
would get a much shallower slope in perfect agreement with 
previous estimates (e.g. Noeske et al. 2007). This result is 
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consistent with Reddy et al. (2012) and cosmological sim¬ 
ulations. Numerical results are consistent with the observa¬ 
tions of Whitaker et al. (2014) and Tomczak et al. (2016) 
for high mass objects (M* > 10'° M©) but there is a severe 
tension at lower masses. As for the higher redshifts consid¬ 
ered in this work, we find that simulations with constant en¬ 
ergy driven winds underpredict the SFR at a fixed mass for 
objects with stellar masses lower than ~ 10'°'^ M© (bottom 
right panel of Figure 5). 

In conclusion, we find that simulations show a good agree¬ 
ment with observations that rely on SED fitting techniques 
for the determination of the intrinsic SFRs and dust correc- 
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tions at z > 1.5. On the contrary, numerical results are not 
consistent with the combination of UV and IR luminosities 
and produce lower star formation rates at a fixed stellar mass 
by almost a factor of 5. This finding confirms the results 
of Katsianis et al. (2015) which addressed higher redshifts 
(z ~ 4 - 7). We demonstrate that the above is true for the 
simulated SFR-M* relations of various groups. It is impor¬ 
tant to note though that the physics assumed for numerical 
modeling are not yet optimized to reflect reality even in the 
state-of-the art simulations and thus it is currently impossi¬ 
ble to determine which observational method produces ro¬ 
bust results. For z < 1.5 numerical results suggest that SFRs 
of high mass objects (log[g(M*/M 0 ) > 10 . 6 ) that rely solely 
on UV luminosities could be underpredicted. This may be 
due to the fact that the dusty environment of high mass 
low redshift galaxies does not allow UV light to escape the 
galaxy and the SFR(/y ends up being underestimated. Ka- 
jisawa et al. ( 2010 ) also suggest that, if a galaxy has star¬ 
forming regions from which one can detect no UV light at 
all due to the heavy dust obscuration (this case occurs more 
frequently for low redshift galaxies with large masses), only 
UV light from relatively less-obscured regions contributes to 
the observed SED, and this results in the underestimation of 
the dust extinction. Furthermore, we see that sample selec¬ 
tion can affect the results for the observed SFR-M* relation. 
The ratio of \o^{S^^i^yman,Ha,hlue Selection!^^^Mass Selection^ tit 
a fixed mass is ~ 0.1 dex to ~ 0.6 dex. The tension increases 
with mass and redshift. The large tension with the normal¬ 
izations and exponents of SFR-M* relations from Lyman- 
break, SFR and Ha selected samples at all redshifts suggests 
that Lyman-break and Ha selections could be biased and 
possibly do not take into account a significant number of ob¬ 
jects. Cosmological simulations predict steeper slopes than 
observations, with an exponent close to unity almost at all 
redshifts. In general, numerical results from different groups 
are in good agreement with mass selected observations. Re¬ 
cent relations that take into account fainter objects are signif¬ 
icantly steeper than those found by past authors and this is in 
agreement with the results of Reddy et al. (2012), Guo et al. 
(2013) and the predictions from cosmological simulations 
(Katsianis et al. 2015). Simulations implemented with vari¬ 
able energy driven and momentum driven winds give similar 
results for the SFR—M* relation and are able to reproduce 
the observables, while models with constant winds fail to 
produce realistic results at low redshift. 

5 Discussion and conclusions 

In the previous sections we have investigated the evolution 
of the SFR-M* relation for galaxies at z ~ 1-4 using cos¬ 
mological hydrodynamic simulations and a compilation of 
observed relations from various groups. In particular, we 
demonstrated that observational studies report a range of 
SFR-M* relations at high redshifts. The reason for the lack 
of consensus is probably related to the fact that different 
groups adopt different methods of selection and calculation 


of the intrinsic properties of galaxies. There has been a con¬ 
siderable effort to constrain the observed star formation rates 
and stellar masses of galaxies but there is not yet a conclu¬ 
sive method of measurements. This has to be addressed by 
future observers since different methods produce different 
results. In addition, in the previous sections we demonstrated 
that the SFR-M* relations obtained from Lyman-break se¬ 
lected and star forming galaxies have typically higher nor¬ 
malizations, when compared with stellar mass based selec¬ 
tion and cosmological hydrodynamic simulations from var¬ 
ious groups. This can be due to the fact that these samples 
are probably biased to include the most luminous-star form¬ 
ing objects and possibly miss a large population of low SFR 
galaxies. This tends to increase the observed mean SFR at 
a fixed mass. The tension between different groups becomes 
smaller at lower redshifts, where observations are more com¬ 
plete. This is possibly due to the fact that various biases re¬ 
lated to sample selection and limits of instrumentation be¬ 
come less severe. Overall, we find that simulations are in 
good agreement with studies which use SED fitting to esti¬ 
mate the intrinsic SERs, dust corrections and stellar masses 
for the observed objects. However, models are unable to re¬ 
produce SER(UVh-IR)-M* relations. 

A comparison between different simulations at z ~ 1 - 4 
suggests that the assumed parameters for the AGN and vari¬ 
able winds feedback implementations do not affect the sim¬ 
ulated SER-M* relation within the range of models con¬ 
sidered. While more extreme models can produce different 
SER-M* relations (Haas et al. 2013), the models in this pa¬ 
per were chosen to reproduce the observed SERE and GSME 
at redshifts z ~ 1 - 7 (Tescari et al. 2014; Katsianis et al. 
2015). We find that the scatter of the SER-M* relation is 
~ 0.2 dex at all redshifts, which is in agreement with esti¬ 
mates from recent observations (Whitaker et al. 2012; Spea- 
gle et al. 2014). 

We note that the implementation of star formation pro¬ 
cesses and ISM physics in simulations has remained essen¬ 
tially unchanged in the last decade. This could be the reason 
of the surprising agreement between older simulations and 
new simulations with an order of magnitude higher resolu¬ 
tions and larger box sizes. Euture numerical codes should 
aim to improve the ISM modeling in order to provide more 
robust estimations of, among other galaxy properties, the 
redshift evolution of the SER-M* relation. 
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A Have past observed SFR-M* relations misguided 
simulations? 

In this appendix we discuss how the uncertainty of dust correction 
laws and selection effects/biases have been affecting the compari¬ 


son with cosmological simulations for the past decade. In figure Al 
we show a compilation of median SFR-M* relations from cosmo¬ 
logical hydrodynamic simulations and semi-analytic models along¬ 
side with observations from various groups. 

Dave (2008) investigated the tension between observed and sim¬ 
ulated results for the SFR-M* and extensively discussed how to 
address it. In general, the numerical results implied steeper rela¬ 
tions with lower normalization than observations. The author con¬ 
sidered various modifications of the theoretical picture of stellar 
mass assembly, but each was found to be in conflict with observa¬ 
tions of high-redshift galaxies. In light of this tension, Dave (2008) 
suggested an evolving IMF to address the discrepancy. Is it pos¬ 
sible though, that the observations used (e.g. Noeske et al. 2007; 
Daddi et al. 2007) to constrain the simulations contained biases 
or/and overestimated the SFR at fixed stellar mass? For redshift 
z ~ 0, we see that the numerical results of Dave (2008) are in ex¬ 
cellent agreement with current state-of-the-art cosmological simu¬ 
lations from the Illustris and EAGLE projects, despite of the fact 
that the resolution is higher and box sizes are larger for the latter 
two simulations. At redshift z ~ 1, the slope of the simulated rela¬ 
tion presented by Dave (2008) is steeper than the observations of 
Noeske et al. (2007). However, in the previous sections we saw in¬ 
dications that the results of Noeske et al. (2007) could have been 
artificially shallower due to the fact that they were not taking into 
acount low mass/SFR objects. The updated results from Guo et al. 
(2013) that take into account low SER objects and 24/jm undetected 
galaxies, are significantly steeper with a power law exponent close 
to unity. Once again, at z ~ 1 we see the perfect agreement be¬ 
tween the numerical results of the Illustris and EAGLE projects and 
the simulated relation from Dave (2008). Moving to z ~ 2, Dave 
(2008) reported a significant tension, with an amplitude offset of 
~ 4 - 5, with the results of Daddi et al. (2007). In the above sec¬ 
tions we saw that the relation given by Daddi et al. (2007) could 
have an artificially high normalization and shallow slope (Bauer 
et al. 2011; Hayward et al. 2014), since the authors preferably se¬ 
lected star forming galaxies and relied on SFRs that were obtained 
from UV-l-IR luminosities. The numerical results of Dave (2008) 
are in good agreement with the simulated SER-M* from the AN¬ 
GUS, Illustris and EAGLE projects, and the mass selected observa¬ 
tions that adopt SED fitting techniques to obtain dust corrections 
and SFRs. Maybe the tension between observed and simulated re¬ 
lations reported by Dave (2008) could have its roots in the fact that 
past observations (Noeske et al. 2007; Daddi et al. 2007) were over¬ 
estimating the SFR at a fixed stellar mass due to methodology, sam¬ 
ple selection effects and biases related to undetected faint objects. 

Sparre et al. (2015) stated that there is good consistency between 
the simulated and observed (Behroozi et al. 2013) SFR-M* rela¬ 
tions at z ~ 0 and z ~ 4. However, at intermediate redshifts there is 
a severe tension. We stress that while the results of Behroozi et al. 
(2013) are a compilation of observations, the authors did not ac¬ 
count for the fact that these observations assumed completely dif¬ 
ferent methods to produce SFR-M* relations. Different methods 
produce completely different results and for this reason they should 
not be compiled directly all together. The results of Behroozi et al. 
(2013) for redshifts z ~ 1 - 2 were mostly based on samples that 
preferably selected star forming galaxies and assumed IR-l-UV lu¬ 
minosities to obtain the intrinsic SFRs (Noeske et al. 2007; Daddi 
et al. 2007; Whitaker et al. 2012). Therefore, the compilation could 
overpredict the SFR at a fixed stellar mass for these redshifts, with 
respect to mass selected surveys that used detailed SED fitting. In 
Eigure Al we show that all models underpredict the SER at a fixed 
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Figure Al. Median values of the SFR-M* relations from different cosmological hydrodynamic simulations for z ~ 0 - 4. The black line is our reference 
model {Ch24.eAJLDW). The dai'k green dotted line is the median fit of the scatter plot presented by Dave (2008). The blue dotted line is the median fit of 
the scatter plot presented by Dutton et al. (2010, SAM). The magenta dashed line is the median line of the scatter plot presented by Sparre et al. (2015, 
lllustris project). The red dashed line is the median line of the scatter plot presented by Furlong et al. (2015, EAGLE project). We cut our SFR(M*) under 
our confidence mass limit of 10^ Mq to make a meaningful comparison with the lllustris and EAGLE projects. There is an excellent agreement between the 
results from cosmological hydrodynamic simulations run by different groups. At each redshift, a panel showing ratios between the different simulations and 
observations with the Ch24.eAjLDW (black solid line) is included. 


stellar mass for z ~ 1-2, with respect to the compilation of obser¬ 
vations in Behroozi et al. (2013). The comparison suggests that it 
is quite possible that the compilation did not take into account the 
faint objects at z 1-2 and therefore suffers from the Malmquist 
bias (the difference with the intrinsic relation suggested by simu¬ 
lations increases at lower masses). A correction to lower star for¬ 
mation rates would be expected, as in Reddy et al. (2012), in order 
to obtain an unbiased relation due to incompleteness. This would 
bring observations and simulations to better agreement. 

In conclusion, the tension reported in the literature likely has its 
roots in the fact that the comparisons have been done using rela¬ 
tions that possibly overpredict the SFR at a fixed stellar mass due 


to selection biases and/or methodology. Cosmological simulations 
in the past decade indicate that the slope of the relation is steeper 
with a lower normalization. Of course, it is possible that this slope 
is the result of a poor representation of physical processes, which 
are implemented in a similar way by simulators (ANGUS, EAGLE, 
lllustris). We demonstrated that the simulated SFR-M* relations 
from various groups are in excellent agreement and largely inde¬ 
pendent of resolution and box size. This somewhat surprising result 
points out how more work is needed to improve the numerical mod¬ 
eling of star formation processes and, in particular, the description 
of the ISM. 
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